Hole-mediated ferromagnetism in tetrahedrally coordinated semiconductors 
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A mean field model of ferromagnetism mediated by delocalized or weakly localized holes in zinc- 
blende and wurzite diluted magnetic semiconductors is presented. The model takes into account: 
(i) strong spin-orbit and kp couplings in the valence band; (ii) the effect of strain upon the hole 
density-of-states, and (iii) the influence of disorder and carrier-carrier interactions, particularly 
near the metal-to-insulator transition. A quantitative comparison between experimental and 
theoretical results for (Ga,Mn)As demonstrates that theory describes the values of the Curie 
temperatures observed in the studied systems as well as explain the directions of the easy 
axis and the magnitudes of the corresponding anisotropy fields as a function of biaxial strain. 
Furthermore, the model reproduces unusual sign, magnitude, and temperature dependence of 
magnetic circular dichroism in the spectral region of the fundamental absorption edge. Chemical 
trends and various suggestions concerning design of novel ferromagnetic semiconductor systems 
are described. 



I. INTRODUCTION 

The^discovery of ferromagnetism in zinc-blendcW'EI III- 
V andtJca II- VI Mn-based compounds allows one to ex- 
plore physics of previously not available combinations 
of quantum structures and magnetism in semiconduc- 
tors. For instance, a possibility of changing the mag=. 
netic phase by light irfl (In,Mn)As/(Al,Ga)Sb andfel 
(Cd,Mn)Te/(Cd,Zn,Mg)Te heterostructures was put into 
the evidence. The injection of spin-polarized carriers 
from (Ga,Mn)As to a (In,Ga)As quantum well in the ab- 
sence of an external magnetic field was demonstrated, 
too.tj It is then important to understand the ferromag- 
netism in these semiconductors, and to ask whether the 
Curie temperatures Tc can be raised to above 300 KJrom 
the present 110 K observed for Gao.947Mno.053As.BtD 

In this paper, we develop theory of the hole-mediated 
ferromagnetism in tetrahedrally coordinated semicon- 
ductors along the lines of a model proposed recently 
by us.El Since we aim at quantitative description of ex-, 
perimental findings, the proposed theoretical approachcl 
makes use of empirical facts and parameters wherever 
possible. We begin the present paper by discussing elec- 
tronic states in p-type magnetic semiconductors. We 
classify the studied systems as charge transfer insula- 
tors, so that our theory is not applicable to materials 
in which d electrons participate in charge transport. We 
note that Mn ions act as both source of localized spins 
and effective mass acceptors. We adapt, therefore, the 
physics of the metal-insulator transition in doped semi- 



conductors for the studied case, and assume that over 
the relevant range of impurity concentrations, the ferro- 
magnetic exchange is mediated by delocalized or weakly 
localized holes. Since the resulting spin-spin coupling is 
long range, we use a mean field approximation to deter- 
mine various thermodynamic, magnetoelastic, and opti- 
cal properties of the system. Particular attention is paid 
to take carefully into account the complex structure of 
the valence band. We then present results of comprehen- 
sive numerical studies which provide qualitative, and in 
many cases quantitative, interpretation of experimental 
finding accumulated over the recent years for (Ga,Mn) As. 
This good agreement between experimental and theoret- 
ical data encourages us to show, in the final part of the 
paper, expected chemical trends, and to propose various 
suggestions concerning the design of novel ferromagnetic 
semiconductor systems. 

In general terms, our results point to the importance 
of the kp and spin-orbit interactions in the physics of the 
hole- mediated ferromagnetism is semiconductors. These 
interactions control the magnitude of Curie temperature, 
saturation value of magnetization, and character of mag- 
netic anisotropies. A comparison of theoretical and ex- 
perimental findings not only emphasizes similarities and 
differences between III-V magnetic semiconductors and 
other ferromagnetic systems, but demonstrates also novel 
aspects of half metallic ferromagnets. The recent compip=. 
hensive reviews present many aspects of IIIfV,El II-VI,E£l 
as well as of IV- VI magnetic semiconductors |L2l which will 
not be discussed here. 
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II. ELECTRONIC STATES IN P-TYPE 
MAGNETIC SEMICONDUCTORS 

A Mn ion in tetrahedrally coordinated 
semiconductors 

We consider zinc-blende or wurzite semiconductor 
compounds, in which the cations are partly substituted 
by magnetic ions, such as Mn. The magnetic ions are 
assumed to be randomly distributed over the cation 
sites, as found by extended x-ray absorption fine strue-. 
ture (EXAFS) studies in the case of Gai-xMn^AstHl 
and Cdi-ajMn^Te.Ej The Mn provides the localized spin 
S = 5/2 and, in the case of III-V semiconductors, acts as 
an acceptor. These Mn acceptors compensate the deep 
antisite donors commonly present in GaAs grown by low- 
temperature molecular beam epitaxy, and produce a p- 
type conduction with metallic resistance for JjheJVLa. con- 
centration x in the range 0.04 < x < 0MiM\!Mi Ac- 
curate analysis of the transport data, complicated by a 
large magnitude of the extraordinary Hall effect, confirms 
the presence of a strong compensation]!!! presumably by 
As antisite donors, as mention above. However, at this 
stage we cannot exclude the existence of self compensa- 
tion mechanisms, such as the formation of Mn AX-like 
centers or donor defects once the Fermi level reaches an 
appropriately deep position in the valence band.tia 

According to optical studies, Mn in GaAs forms an 
acceptor, center characterized by a moderate binding 
energyE3 E a = 110 meV, and a small magnitude of the 
energy differenGe^between the triplet and singlet state of 
the bound hole£3E3 Ae = 8 ± 3 meV. This small value 
demonstrates that the hole introduced by the divalent 
Mn in GaAs does not reside on the d shell or forms 
a Zhang- Rice-likej-singlctJ^'Ej but occupies an effective 
mass Bohr orbit BEj Thus, due to a large intra-site corre- 
lation energy U, (Ga,Mn)As can be classified as a charge- 
transfer insulator, a conclusion consistent with photoe- 
mission spectroscopy.EiiEj At the same time, the p-d hy- 
bridization results in a spin-dependent coupling between 
the holes and the Mn ions, H p d = —(3N sS. Here /3 is the 
p-d exchange integral and N Q is the concentration of-thn 
cation sites. The-analysis of both photoemission datacj'E£l 
and magnitude^ of Ae leads to the exchange energy 
f3N ps —1 eV. Similar values of (3N a are observed in 
II- VI diluted majpj|etic semiconductors with comparable 
lattice constants.E3 This confirms Harrison's suggestion 
that the hybridization matrix-elements depend primarily 
on the inter-atomic distanceEO According to the model in 
question, the magnetic electrons remain localized at the 
magnetic ion, so that they do not contribute to charge 
transport. This precludes Zener's double exchangeEH as 
the mechanism leading to ferromagnetic correlation be- 
tween the distant Mn spins. At the same time, for some 
combinations of transition metals and hosts, the "chemi- 
cal" and exchange attractive potential introduced by the 
magnetic ion can be strong enough to bind the hole on 



a local orbitJ2M2l In an intermediate regime, the proba- 
bility of finding the hole around the magnetic ion is en- 
hanced, which results in the apparent increase of |/3iV | 
with decreasing xEj 

B Two-fluid model of electronic states near the 
metal-insulator transition 

Ionized impurity and magnetic scatterings lead to lo- 
calization of the effective mass holes introduced by Mn 
in III-V compounds or by acceptors in the case of II- 
VI materials. It is, therefore, important to discuss the 
effect of Anderson-Mott localizatiop-|On the onset of fer- 
romagnetism. The two-fluid modeO constitutes the es- 
tablished description of electronic states in the vicinity 
of the Anderson-Mott metal-insulator transition (MIT) 
in doped semiconductors. According to that model, the 
conversion of itinerant electrons into singly occupied im- 
purity states with increasing disorder occurs gradually, 
and begins already on the metal side of the MIT. This 
leads to a disorder-driven static phase separation into two 
types of regions: one populated by electrons in extended 
states, and another that is totally depleted from the elec- 
trons or contains singly occupied impurity-like states. 
The latter controls the magnetic response of doped non- 
magnetic semiconductorsE] and gives rise to the pres- 
ence of bound magnetic polarons (BMRljmi-both sides 
of the MIT in magnetic semiconductors.c3£3'U Actually, 
the formation of BMP .shifts the MIT towards the higher 
carrier concentrationsrSEH On crossing the MIT, the 
extended states become localized. However, according to 
the scaling theory of the MIT, their localization radius 
£ decreases rather gradually from infinity at the MIT 
towards the Bohr radius deep in the insulator phase, so 
that on a length scale smaller than £ the wave function re- 
tains an extended character. Such weakly localized states 
are thought to determine the static longitudinal and Hall 
conductivities of doped semiconductors. The central sug- 
gestion of the recent modeltJ is that the holes in the ex- 
tended or weakly localized states mediate the long-range 
interactions between the localized spins on both sides of 
the MIT in the III-V and II- VI magnetic semiconductors. 

As will be discussed below, the Curie temperature 
Tq is proportional to the thermodynamic spin density- 
of-states p s which, in turn, is proportional to the spin 
susceptibility of the carrier liquid \s- Like other ther- 
modynamic quantities, p s does not exhibit any critical 
behavior at the MIT. However, p s exhibits large space 
fluctuations at criticality, which will result in local fluc- 
tuations of magnetic properties. The quantitative renor- 
malization of p s by disorder will depend on its micro- 
scopic nature, for instance, on the degree of compensa- 
tion. The enhancement of p s by the carrier-carrier inter- 
actions can be described by the Fermi-liquid parameter 
Ap, p s — > App s £3 The value of Ap = 1.2, as evaluated^ 
by the local-spin-dcnsity approximation for the relevant 
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hole concentrations, has been adopted for our compu- 
tations. We note that disorder-modified carrier-carrier 
interactions in the triplet particle-hole channel tend to 
enhance Ap, which maypeven drive the system towards 
a Stoner-like instabilityES In magnetic systems, how- 
ever, spin-flip scattering by fluctuations of magnetization 
makes this enhancement mechanism rather inefficient. 

The two-fluid model is consistent with the recent EPR 
results) 3 ^ which point the coexistence of the neutral Mn 
acceptors and ionized Mn d 5 states in some f$fize of Mn 
concentration. Furthermore, the observeaOH sign of 
magnetic circular dichroism in (Ga,Mn)As suggests the 
presence of the Fermi-liquid like states on th e both sides 
of the MIT, as we shall discuss in Sec. IV D| . 
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C Valence band structure and exchange splitting 
of the hole subbands 

Since the valence band originates merely from the an- 
ion p and cation d wave functions, the exchange interac- 
tion mediated by the holes is expected to be strongly 
affected by anisotropy of hole dynamics and the cou- 
pling between the spin and orbital orbit degrees of free- 
dom. To take those effects into account, the hole dis- 
persion and wave functions are computed by diagonaliz- 
ing the 6x6 Kohn-Luttinger kp matrixHj In this model, 
four Tg and two T7 bands are taken explicitly into ac- 
count, whereas other bands are included by second order 
perturbation theory. The model is developed for zinc- 
blende and wurzite semiconductors. It allows for warp- 
ing, quantizing magnetic fields, and biaxial strain but no 
terms associated with the lack of the symmetry inversion 
are taken into account. The effect of the spin-dependent 
interaction between the holes and the Mn spins is de- 
scribed in terms,©f the virtual-crystal and molecular-field 
approximations,^! so that 



H pd = l3sM(r)/gn B , 



(1) 



where M(r) is the magnetization of the localized spins 
that carry magnetic moment —SgfiB, where S — 5/2 and 
.9-2.0. 

The explicit form of the kp 6x6 matrices, together 
with the matrix H p d derived by us in the Luttinger- 
Kohn representation for arbitrary directions of magne- 
tization M, are displayed in Appendix A. We note that 
for the case under consideration, involving large exchange 
interaction and high kinetic energy, particularly impor- 
tant are off diagonal terms describing the p-d coupling 
between the Tg and bands. A numerical procedure 
that serves to determine the concentration, free energy, 
wave functions, and optical characteristics of the holes 
is outlined in Appendix B. The adopted values of the 
Luttinger parameters ji and the spin-orbit splittings A G 
are summarized in Appendix C for various III-V and 
II- VI parent compounds. According to photoemission 
studiesEJ (3N Q = -1.2 ± 0.2 eV for (Ga,Mn)As. The 



FIG. 1: Fermi energy as a function of the hole concen- 
tration p computed from the 6x6 Luttinger model for 
kp parameters displayed in the figure. Inset shows cross 
section of the Fermi sphere for p — 3.5 x I0 20 cm -3 . 

value f3N = —1.2 eV is taken for the computations, 
though future works may reveal some dependence of [3 of 
Eq. (nl) on Mn and/or hole concentrations because of en- 
ergetic proximity of the Fermi and relevant Mn d-levels. 
Such a dependence could result also from corrections to 
the virtual-crystal and molecular-field approximations if 
a short-range part of the Mn potential is attractive for 
the valence band hole.E3 

Figure [l| presents the dependence of the low tempera- 
ture Fermi energy ep on the hole concentration p com- 
puted for the band parameters corresponding to GaAs. 
No effects of p-d exchange, disorder, or Coulomb inter- 
actions between the holes are taken into account. The 
corresponding density-of-states (DOS) effective mass in- 
creases from 0.67m o in the limit of small hole concentra- 
tions to the value of 1.35m for p — 5 x 10 20 cm~ 3 . This 
increase is caused by the kp interaction between the Tg, 
and T7 bands, which-because of a relatively small mag- 
nitude of A in GaAs-is important for the relevant hole 
densities. The inset to Fig. [j] shows the cross section of 
the Fermi surface for p — 3.5 x 10 20 cm -3 , which corre- 
sponds to sf — —195 meV in respect to energy of the T s 
point. Two valleys, the heavy and light hole subbands, 
are visible. 

A strong and complex influence of the p-d interaction 
and strain upon the valence band is shown in Fig. ^|. The 
cross sections of the Fermi spheres are depicted for p — 
3.5 x I0 20 cm' 3 , which now corresponds to ep ~ —165 
meV, and for the parameter of the exchange splitting, 

B G = A F pM/6gfx B , (2) 

taken as B G = -30 meV. If (3N = -1.2 eV and 
Ap = 1.2, this magnitude of Bq occurs for the saturation 
value of magnetization M at x — 0.05. We note at this 
point that some of the effects, e.g., the direction of spin- 
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polarization vector depends on the sign of (3 and thus Bq, 
whereas others, like the Curie temperature, are propor- 
tional to P 2 . In the absence of strain, e = 0, the four fold 
degeneracy at the T point is lifted by the p-d exchange, 
and the corresponding energies at k = are ±3£?g and 
±Bq for A D 3> \Bq\- However, the splitting at non-zero 
wave vectors depends on the relative orientation of M 
and fco In particular, since the spin of the heavy hole is 
polarized along k for e = 0, the exchange splitting is seen 
to vanish for k _L M. This mixing of orbital and spin 
degrees of freedom, together with highly nonparabolic, 
anisotropic, and mutually crossing dispersion relations 
constitute important aspects of the hole-mediated ferro- 
magnetism in terahcdrally coordinated semiconductors. 



III. HOLE-INDUCED FERROMAGNETISM IN 
SEMICONDUCTORS 

A Short-range atiferromagnetic superexchange and 
ferromagnetic double exchange 

In addition to the interaction between the carriers 
and localized spins, the p-d hybridization leads also 
to the superexchange, a short-range antiferromagnetic 
coupling between the Mn spins. The superexchange 
is mediated by spin-polarization of occupied electron 
bands, in contrast to the Zener ferromagnetic exchange 
which is mediated by spin-polarization of the carrier liq- 
uid. The antiferromagnetiCj-exchange dominates in un- 
doped II- VI semiconductors,!^] and also in a compensated 
Gai-xMn^AsiSn.EJ In order to take the influence of this 
interaction into account, it is convenient to parameterize 
the dependence of magnetization on the magnetic field 
in the absence of the carriers, M a (H), by the Brillouin 
function B5 according to 

M (H) = gn B SN x ef fBs[g»BH/k B (T + T AF )}, (3) 

where two empirical parameters, the effective spin con- 
centration x e ffN a < xN a and temperature T e ff > T, 
take the-uresence of the superexchange interactions into 
account Ejo The dependencies x e ff(x) and Taf(%) 
7, . f(x) — T are knownErc2l for II- VI compounds. 

In order to elucidate the effect of doping on x e ff and 



Taf we refer to the two-fluid model described in Sec. II B 



In terms of that model the delocalized or weakly local- 
ized holes account for the ferromagnetism. Actually, the 
participation of the same set of holes in both charge 
transport and-the ferromagnetic interactions is shown, 
in (Ga,Mn)AfJiil and in (Zn,Mn)TeJiS by the agreement 
between the temperature and field dependencies of the 
magnetization deduced from the extraordinary Hall ef- 
fect, Mjj, and from direct magnetization measurements, 
Mjj, particularly in the vicinity of Tq- However, below 
Tc and in the magnetic fields greater than the coercive 
force, while Mh saturates (as in standard ferroisaaenets), 
Mb continues to rise with the magnetic fieldJi-fi Since 



Mh is proportional to spin polarization of the carriers, its 
saturation may reflect a saturation of hole polarization, 
which-for appropriately low values of the Fermi energy- 
can occur even if the Mn spins are not totally spin polar- 
ized (half metallic case). It appears, however, that Mn 
spins in the regions depleted from the carriers (which, 
therefore, do not participate in the long range magnetic 
order) contribute also to a slowly increasing component 
ofMo(ff). 

According to the two fluid model, part of the carriers 
are trapped on strongly localized impurity states, and 
thus form bound magnetic polarons (BMP). If there is 
an exchange coupling between the two fluids, the BMP 
participate in the formation of the ferromagnetic order. 
Furthermore, the coupling between the BMP appears 
to be ferromagnetic, at least in some range of relevant 
parametersOc3 To gain the Coulomb energy, the BMP 
are preferentially formed around close pairs of ionized 
acceptors. In the case c±4II-V materials this leads, via 
Zener's double exchange,!^ to a local ferromagnetic align- 
ment of neighbor Mn d 5 negative ions,c3so that x « x e ff 
and Taf ~ 0. By contrast, in II- VI compounds, for which 
acceptor cores do not carry any spin and the degree of 
compensation is low, BMP are not preferentially formed 
around Mn pairs, so that the close pairs remained anti- 
ferromagnetically aligned, even in p-type samples. The 
presence of a competition between the ferromagnetic and 
antiferromagnetic interactions in II- VI compounds, and 
its absence in III-V materials, constitutes the important 
difference between those two families of magnetic semi- 
conductors. 



B Zener model of ferromagnetic interactions 
mediated by free carriers 

As mentioned above, we assume that weakly localized 
or delocalized holes mediate long-range, ferromagnetic in- 
teraction between the spins. ZenereZI first proposed a 
model of ferromagnetism driven by the exchange cou- 
pling of the carriers and the localized spins. According 
to that model, spin-polarization of the localized spins 
leads to spin-splitting of the bands, which results in the 
lowering of the carrier energy. At sufficiently low tem- 
perature, this lowering overcompensates the increase of 
the free energy caused by a decrease of entropy, which 
is associated with the polarization of the localized spins. 
However, the Zener model was later abandoned, as nei- 
ther the itinerant character of the magnetic electrons nor 
the quantum (Friedel) oscillations of the electron spin- 
polarization around the localized spins were taken into 
account, both of these are now established to be critical 
ingredients for the theory of magnetic metals. In partic- 
ular, a resulting competition between ferromagnetic and 
antiferromagnetic interactions in metals leads rather to 
a spin-glass than to a ferromagnetic ground state. In 
the case of semiconductors, however, the mean distance 
between the carriers is usually much greater than that 
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between the spins. Under such conditions, the exchange 
interaction mediated by the carriers is ferromagnetic for 
most of the spin pairs, which reduces the tendency to- 
wards spin-glass freezing. Actually, for a random distri- 
bution of the localized spins, the mean-field value of the 
Curie temperature Tc deduced from the Zener model is 
equal to that obtained from the RudeimarL_Kittel, Ka- 
suya, and Yosida (RKKY) approach, l3c3'EiE3 in which 
the presence of the Friedel oscillations is explicitly taken 
into account. 




c) 



k [10 7 cm"'] 



C Mean-field model of Curie temperature and 
thermodynamic properties 

As reported elsewherej^ the Zener model describes cor- 
rectly the experimental values of Tc in both (Ga,Mn)As 
and (Zn,Mn)Tc, provided that band structure effects are 
taken into account. The starting point of the model is 
the determination how the Ginzburg-Landau free-energy 
functional F depends on the magnetization M of the lo- 
calized spins. The hole contribution to F, F C [M] is com- 
puted by diagonalizing the 6x6 Kohn-Luttinger matrix 
together with the p-d exchange contribution, and by the 
subsequent computation of the partition function Z, as 
described in Appendix B. This model takes the effects 
of the spin-orbit interaction into account, a task diffi- 
cult technically within the RKKY approach, as the spin- 
orbit coupling leads to non-scalar terms in the spin-spin 
Hamiltonian. Moreover, the indirect exchange associated 
with the virtual spin excitations between the valence sub- 
bands, the Bloembergen-Rowland mechanisms] 2 ^ is auto- 
matically included. The remaining part of the free energy 
functional, that of the localized spins, is given by 



F S [M] = / dM H(M ), 
Jo 



(4) 



where H(M ) is given in Eq. (||). By minimizing F[M] = 
F C [M] + F S [M] with respect to M at given T, H, and 
hole concentration p, one obtains M(T, H) as a solution 
of the mean-field equation, 



M = Li g iJ. B SN x e ffB s 



gH B (-dF c [M]/dM + H) 



k B (T + T AF ) 



(5) 

Near the Curie temperature Tc, where M is small, we 
expect F C [M] — F c [0] ~ M . It is convenient to parame- 
terize this dependence by a spin density-of-states p s , 



F C [M] = F c [0] - A F p s p 2 M 2 /2(2gp B ) 



(6) 



The spin density-of-states p s is related to carrier mag- 
netic susceptibility according to x = A F {g* p, B ) 2 p s / A 
and, in general, has to be determined numerically by 
computing i^ c [M]. By expanding Bg(Af) we arrive to, 

T c = x eff N S(S +l)P 2 A F p s (T c )/Uk B - T AF . (7) 



FIG. 2: Cross sections of the Fermi spheres for saturation 
value of magnetization M in GaogsMno.osAs with hole 
concentration p = 3.5 x 10 20 cm -3 . Results for various 
orientations of magnetization M and biaxial strains e xx 
are shown: (a) M |j [100]; e xx = 0. (b) M || [100]; e xx = 
-2%. (c)M || [001]; e xx = +2%. A strong dependence 
of the splitting on the relative orientation of the wave 
vector and magnetization is visible. 



For a strongly degenerate carrier liquid — e F \/k B T ^S> 1, 
as well as neglecting the spin-orbit interaction, p s be- 
comes equal to the total density-of-states p for intra-band 

2fc2 



charge excitations, where p = m* DOS k F /it h 

In order to check the quantitative significance of car- 
rier entropy, the computed values of Tc were compared 
to those obtained assuming strong degeneration of the 
carrier liquid. As shown in Fig. |3j such an assumption 
leads to error smaller than 1% if e F \/k B T > 10. Thus, 
in this range, the carrier energy, not free energy, was used 
for the evaluation of Tc- Furthermore, we take £ e // = x 
and Taf = 0, as discussed above. 

The values of Tc as a function of the hole concentra- 
tion p for Gao.95Mno.05As, as computed by our model, 
are shown by the solid line in Fig. ^. Since the magni- 
tude of Tc is directly proportional to x e ffN (3 2 A F , the 
theoretical results can be easily extended to other values 
of x, (3 or A F . In particular, Tc — 300 K is expected for 
x = 0.125 and p = 3.5 x 10 20 cm" 3 . 



D Validity of the model and comparison to other 
approaches 

Equation (0) with p s — p, that is neglecting effects of 
the spin-orbit interaction, has already. Jaeen derived by 
a number of equivalent methods oaoa In the present 
work F c [Af], which served to determine Tc, M(T,H), 
anisotropy fields, and optical spectra was obtained from 
the 6x6 kp model, as described in Appendices A and B. 
In order to illustrate the importance of band structure 
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FIG. 3: The computed ratio the actual Curie tempera- 
ture Tc to that obtained assuming that the hole liquid 
is strongly degenerate as a function of the degeneration 
factor. 



effects, the evaluation of Tc has been performed employ- 
ing various models. The two straight dashed lines in 
Fig. U depict the expected Tc(p) if the warping is ne- 
glected (72 =73) as well as a large and small spin-orbit 
splitting Ep s — Ep 7 = A Q is assumed, respectively. We 
see that the above assumptions are approximately ful- 
filled in low and high concentration range, respectively, 
but the calculation within the full model is necessary in 
the experimentally relevant region of iJpje .intermediate 
hole densities. In particular, 4x4 modeOcj (A D ^ \£f\) 
would be reasonably correct for (Zn,Mn)Te, where A D is 
about 1 eV but not for (Ga,Mn)As, for which A D is al- 
most three times smaller. Finally, the dotted line shows 
Tc evaluated replacing p s by p. We see that p ~ p s al- 
ready for p 10 20 cm" 3 , despite that \ef\ is still two 
times smaller than A„. In contrast, p s is more than two 
times smaller than p in the limit of small p. This reduc- 
tion of the spin susceptibility, and thus of Tc, stems from 
the absence of the exchange splitting of the heavy hole 
subband for k _L M. 

As described above, Tc can be computed by minimiz- 
ing the free energy, and without referring to the explicit 
form of the Kohn-Luttinger amplitudes Uik- Since near 
Tc the relevant magnetization M is small, the carrier free 
energy, and thus Tc, can also be determined from the lin- 
ear response theory. The corresponding p s assumes the 
standard form 
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FIG. 4: Curie temperature as a function of the hole con- 
centration for Gao.95Mno.05As computed from the 6x6 
Luuttinger model (solid line). Straight dashed lines rep- 
resent results obtained assuming a large and small value 
of the spin-orbit splitting A OJ respectively. The dotted 
line is calculated neglecting the effect of the spin-orbit 
interaction on the hole spin susceptibily. 
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\(ui,k\sM\u 3} k+ q )\ 2 fi(k)[l - fjjk + q)] 
E. ] {k + q)-E l {k) 



(8) 

where sm is the component of spin operator along the 
direction of magnetization and fi(k) is the Fermi-Dirac 
distribution function for the i-th valence band subband. 
The equivalence of p s as given by Eqs. (g) and (||) can be 
checkad for the 4x4 spherical model by using the explicit 
formed of Ei(k) and Such a comparison demonstrates 
that almost a 30% contribution to Tc originates from 
interband polarization. 

It is straightforward to generalize the model for 
the case, uf .the carriers confined to the d-dimensional 
space-LjOo The tendency towards the formatioa of 
spin-density waves in low-dimensional systemsE3E3 as 
well as possible spatial correlation in the distribution of 
the magnetic ions can also be taken into account. The 
mean- field value of the critical temperature T q , at which 
the system undergoes the transition to a spatially mod- 
ulated state characterized by the wave vector q, is given 
by the solution of the equation, 

P 2 A F (q,T q )p s (q,T q ) J dCXo(q, T q , C)\MC)\ 4 = 

Here q spans the d-dimensional space, <fi (C) is the en- 
velope function of the carriers confined by a (3 — d)- 
dimensional potential well V(£); g and \o denote the 
Lande factor and the q-dependent magnetic susceptibil- 
ity of the magnetic ions in the absence of the carriers, re- 
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spectively. Within the mean-field approximation (MFA) , 
such magnetization shape and direction will occur in the 
ordered phase, for which the corresponding T q attains 
the highest value. A ferromagnetic order is expected in 
the three dimensional (3D) case, for which a maximum 
of p s (q) occurs at q = 0. 

Since within the MFA the presence of magnetization 
fluctuations is neglected, our model may lead to an over- 
estimation of the magnitude of M, and thus of Tq- 
To address this question we recall that the decay of 
the strength of the carrier- mediated exchange interaction 
with the distance between the spins r is described by the 
RKK¥,|f-U|iaction, which in the 3D situation assumes the 
formula 

J(r) - [sin(2/fc F r) - 2k F r cos(2fc F r)]/(2fc F r) 4 . (10) 

In the case of semiconductors an average distance be- 
tween the carriers r c = (Anp/3)^ 1 ^ 3 is usually much 
greater than that between the spins r$ = (47raiV /3) _1 / 3 . 
This means that the carrier-mediated interaction is fer- 
romagnetic and effectively long range for most of the 
spins as the first zero of the RKKY function occurs at 
r w 1.17r c . A theoretical studyEil of critical exponents 
for a (i-dimensional space showed that as long as a < d/2 
in the dependence J(r) ~ l/r d+a , the mean-field ap- 
proach to the long wavelength susceptibility x(T) is valid, 
a conclusion not affected presumably by disorder in the 
spin distribution. At the same time, both relevant length 
scale r c , not rs, and the critical exponents^ rj = 2 — a 
and v = I/it, point to much faster decay of xil) with q 
than that expected from the classical Ornstein-Zernike 
theory.LJ This indicates that the MFA should remain 
valid down to at least \T - T c \/T c » rs/r c <C 1. Actu- 
ally, the decay of x(<l) with q in the range < q < 2Uf is 
corroborated by the observation of smaller critical scat- 
tering of the holes by magnetisation fluctuations than 
that calculated for x{l) — X(0).t3 

Recently, Monte-Carlo studies of carrier-mediated- 
romagnetism in semiconducors have been initiated. 
Such an approach has a potential to test the accuracy of 
the MFA and to determine the actual spin configuration 
corresponding to the ground state. Prelipiuiary results 
appear to confirm the validity of the MFAEjIiH as well as 
to indicate a possibility of the existence of non-caLLinear 
magnetic structures in low-dimensiopal systems.LZI An- 
other significant recent developments is the examina- 
tion of spin-wave excitations, their spectrum and ef- 
fect on magnetization. A strong reduction of Tc was 
predicted,E3 though it should be noted that the spin wave 
approximation breaks usually down at criticality. In con- 
trast, such an approach offers valid results at low tem- 
peratures, provided that effects of magnetic anisotropy 
are thoroughly taken into account. Furthermore, the 
existence of ferromagnetic ground state was confirmed 
by ab initio total energy computations for magnetic III- 
V compounds E3H'E2l However, to what extend the em- 
ployed procedures have been capable to capture accu- 



rately correlation effects on the open d-shells seems to be 
unclear by now. 

Finally, we would like to stress once more that if con- 
centrations of the carriers and the spins become compara- 
ble, r c < rs, randomness associated with the competition 
of ferromagnetic and antiferromagnetic interactions can 
drive the system towards a spin-glass phaseH In the case 
of II- VI compounds, the antiferromagnetic component is. 
additionally enlarged by the superexchange interactions 
Furthermore, scattering by ionized impurities, and the 
associated non-uniform distribution of carriers in semi- 
conductors near the MIT, may enhance disorder effects 
even further. We note also that in the extreme case, 
r c *C rs, the Kondo effect, that is dynamic screening of 
the localized spins by the sea of the carriers may preclude 
both ferromagnetic and spin-glass magnetic ordering. 



IV. COMPARISON OF THEORETICAL AND 
EXPERIMENTAL RESULTS FOR (GA,MN)AS 

A Curie temperature and spontaneous 
magnetization 

The most interesting property of Gai-^Mn^As epilay- 
ers is the large magnitudaxif Tc, up to 110 K for the Mn 
concentration x = 5.3%£ll3 Because of this high Tc, the 
spin-dependent extraordinary contribution to the Hall re- 
sistance Rh persists up to 300 K, making, an .accurate 
determination of the hole density difficult .Qtialij'Ej How- 
ever, the recent measurementE 2 ! of Rh up to 27 T and at 
50 mK yielded an unambiguous value of p — 3.5 x 10 20 
cm~ 3 for the metallic Gag. 947Mno.o53As sample, in which 
Tc = 110 K is observedE The above value of p is about 
three times smaller than xN Q , confirming the importance 
of compensation in (Ga,Mn)As. As shown in Fig. [|, the 
numerical results lead to Tc — 120 K for x = 0.05, and 
thus, T c = 128 K for x = 0.053 and p = 3.5 x 10 20 
cm~ 3 . We conclude that the proposed model describes, 
with no adjustable parameters, high values of Tc found 
in (Ga,Mn)As. Furthermore, scaling t heor y of electronic 
states near the MIT, discussed in Sec. [IB, makes it pos- 



sible to explain the presence of the ferromagnetism on the 
both sides of the MIT, and a non-critical evolution of Tc 
across the critical point u A comparison between theoret- 
ical and experimental data in a wider range of Mn and 
hole concentrations requires reliable information on the 
hole density in particular samples, which is not presently 
available in the case of (Ga,Mn)As. 

In the case of (Zn,Mn)Te:N, the hole concentration can 
readily be, determined by the Hall effect measurements 
at 300 K.LTCJ The absence of the extraordinary Hall ef- 
fect at 300 K stems form the much lower values of Tc, 
0.75 < T c < 2.4 K for 5% > x > 2% and 10 19 < p < lO 2 " 
cm~ 3 . The model in question explains satisfactorily^ 
Tc(x,p). Three effects conspire to make Tc much greater 
in p-(Ga,Mn)As than in p-(Zn,Mn)Te at given p and x. 
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First, as already mentioned, the small value of the split- 
ting A Q makes the reduction of Tq by spin-orbit coupling 
of minor importance in the case of (Ga,Mn)As. Sec- 
ond, because of a smaller lattice parameter, the prod- 
uct (3N a is greater in (Ga,Mn)As. Finally, ferromagnetic 
double-exchange between closely lying pairs of Mn ions 
is stronger than antiferromagnetic superexchange in com- 
pensated (Ga,Mn)As. In contrast, antiferromagnetic su- 
perexchange remains significant in p-(Zn,Mn)Te, as the 
Mn atoms are electrically inactive in II- VI compounds. 

Important characteristics of ferromagnets are the mag- 
nitude and temperature dependence of the spontaneous 
magnetization M below Tc- Experimental studies of 
metallic samples indicates that on lowering temperature 
M increases as the Brillouin function reaching the sat- 
uration value M s for T — > 0.0 This behavior indicates 
that a molecular field H* acting on the Mn spins is pro- 
portional to their magnetization M. Figure ^| presents 
values of M/M s as a function of T/Tc calculated with 
no adjustable parameter for (Ga,Mn)As containing vari- 
ous hole and Mn concentrations. We see that indeed the 
computed Mn spin magnetization M(T)/M S follows, to 
a good approximation, the Brillouin function, except for 
materials with rather small values p or large values of 
x. The latter cases correspond to the half-metallic sit- 
uation, for which only the ground-state hole subband is 
populated even at M < M s , so that the molecular field 
produced by the carriers H* attains its maximum value 
and, therefore, ceases to be proportional to M . Never- 
theless, the saturation value M s can be reached provided 
that temperature is sufficiently low, ksTeff <C Sg^sH* . 
Such a half-metallic behavior is observed in the case of 
ferromagnetic correlation imppaed by a dilute hole liquid 
in (Cd,Mn)Te quantum wellsE3 

B Hole magnetization and spin polarization 

The magnitude of magnetization presented in the pre- 
vious section was computed neglecting a possible contri- 
bution originating from hole magnetic moments. Such a 
contribution can be significant as the hole liquid is spin- 
polarized for non-zero magnetization of the Mn spins. 
Because of the spin-orbit interaction, the hole magneti- 
zation consists of two components. One comes from mag- 
netic moments of the hole spins, described by the Lande 
factor of the free electrons and the Luttinger parameter k. 
Another contribution, absent for localized carriers, orig- 
inates from diamagnetic currents, whose magnetic mo- 
ments can be oriented along the spin polarization by the 
spin-orbit interaction. The evaluation of the latter re- 
quires the inclusion of the Landau quantization in the kp 
hamiltonian. The carrier magnetization M c is then given 
by, 
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FIG. 5: The computed emerging of spontaneous mag- 
netization M below the Curie temperature Tc for vari- 
ous hole concentrations p (a) and Mn contents x (b)in 
Gao.95Mno.05As (points). If the exchange splitting of 
the valence band is much smaller than the Fermi energy 
(large p, small x), the evolution of magnetization follows 
the Brillouin function (dashed lines). 



M c = - lim dG c 



(11) 
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where the Gibbs thermodynamic potential is calculated 
for a given value of the Mn spin magnetization M and 
the Fermi energy sf(M) as a function of the magnetic 
field acting on the carriers H. In general, the eigenvalue 
problem for the holes in the magnetic field cannot be 
transformed into an algebraic equation. Such a transfor- 
mation is, however, possible if the warping is neglected. 
We have, therefore, calculated the hole magnetization M c 
disregarding the anisotropy, that is assuming 71 = 6.85, 
72 = 73 = 2.58, and k = 1.2. The explicit form of the 
relevant Luttinger matrices is displayed in Appendix A. 
The partition function Z was computed by summing over 
the Landau index n, the wave vector k z , and the six hole 
subbands. 

The results of the computations are shown in 
Fig. ^, where M c is plotted as a function of p for 
Bq — Ajrf3M/6gfiB = —30 meV. In this paper, 
for sake of comparison with experimental results we 
depict magnetization in the SI units according to 
/j, o M[T]=47rl0 _4 M[emu]. It is seen the diamagnetic (or- 
bital) contribution to M c is negative and dominant. The 
spin term is positive, which for the antiferromagnetic sign 
of the p-d exchange integral f3 points to the negative sign 
of the hole Lande factor gh . A visible decrease of the spin 
contribution for the large p corresponds to a cross-over 
to the free electron value g — 2.0 occurring when \ep\ 
approaches A Q . In Fig. |7|, M c is plotted versus M for 
various p. It is seen that for the employed values of the 
parameters, M c reaches only 5% of /j. M = 65 mT, which 
corresponds to the saturation value of Mn magnetization 
for x = 0.05. A rather weak magnitude of M c results 
from a partial compensation of the spin and orbital con- 
tributions to M c as well as from smaller concentration 
and spin of the holes in comparison to those of the Mn 
ions. We conclude that delocalized or weakly localized 
holes give a minor contribution to the total magnetiza- 
tion. Accordingly M c is neglected when determining the 
direction of easy axes and the magnitude of anisotropy 
fields. P 

In view of the on-going experimentscl on electrical spin 
injection from (Ga,Mn)As, an important question arises 
what is the degree of hole spin polarization V as a func- 
tion of p and Bq. Furthermore, V appears to control the 
magnitude of the extraordinary Hall effect. It is, there- 
fore, interesting to determine conditions, under which the 
usual assumption about the linear relation between V 
and magnetization of the Mn spins M is fulfilled. 

The contribution of all four hole subbands to the Fermi 
cross section visible in Fig. || indicates that the exchange 
splitting is too small to lead to the total spin polariza- 
tion for p = 3.5 x 10 20 cm~ 3 and x = 0.05 (\B G \ = 30 
meV). Furthermore, a destructive effect of the spin-orbit 
interaction on the magnitude of V can be expected. In 
order to evaluate V = 2(sm)/p we note that according 
toEq. d), 
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FIG. 6: Magnetization M c of the hole liquid (squares) in 
Gai-^Mn^As computed as a function of the hole concen- 
tration for the spin splitting parameter Bq = — 30 meV 
(the latter corresponds to the saturation value of Mn spin 
magnetization M for x = 0.05). Crosses show spin and 
orbital contribution to M r . 
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FIG. 7: Computed hole magnetization M c as a function 
of Mn spin magnetization M for various hole concentra- 
tions in Ga 1 _ 2 ,Mn a .As. 



V = 



2 gf i B dF c (M) 
(3p dM 



(12) 
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SPLITTING PARAMETER, \BJ [meV] 

FIG. 8: Computed degree of spin polarization of the hole 
liquid as a function of the spin splitting parameter for 
various hole concentrations in Gai-^Mn^As (Bq = —30 
meV corresponds to the saturation value of Mn spin mag- 
netization for x = 0.05). The polarization of the hole 
spins is oriented in the opposite direction to the polar- 
ization of the Mn spins. 

Figure || presents the dependence of V on Bq for the 
experimentally relevant range of p. We see that \T\ 
tends to saturate with \Bq\, and thus with M. This 
means that for large values of the splitting Bq, magne- 
tization M evaluated from the extraordinary Hall effect 
is underestimated. At the same time, the calculation 
demonstrates that despite the spin-orbit interaction \V\ 
becomes greater than 0.8 for 3|£?g| > \sf\- This is due to 
the fact that the redistribution of the holes over the spin 
subbands occurs in the way which maximizes the gain 
of the exchange energy, and thus the magnitude of \V\. 
This in contrast to the case \ef\ S> \Bq\, for which the 
spin polarization is reduced by a factor greater than two 
from the value corresponding to the absence of the spin- 
orbit coupling at low hole concentrations. We also note 
that because of the antiferromagnetic character of the p- 
d coupling ((3N < 0) , the polarization of the hole spins 
is oriented in the opposite direction than the polarization 
of the Mn spins. 

C Easy axis and anisotropy field 

Already early studies of a ferromagnetic phase in 
(Ga,Mn)As epilayers demonstrated the existence of sub- 
stantial magnetic anisotropyEa Magnetic anisotropy is 
usually associated with the interaction between spin and 
orbital degrees of freedom of the magnetic electrons. Ac- 
cording to the model in question, these electrons are in 
the d 5 configuration. For such a case the orbital mo- 
mentum L — 0, so that no effects stemming from the 
spin-orbit coupling are to be expected. To reconcile the 



model and the experimental observations, we note that 
the interaction between the localized spins is mediated by 
the holes, characterized by a non-zero orbital momentum. 
An important aspect of the present model is that it does 
take into account the anisotropy of the carrier-mediated 
exchange interaction associated with the spin-orbit cou- 
pling in the host material, an effect difficult to include 
within the standard approach to the RKKY interaction. 

We start by considering an unstrained thin film with 
the [001] crystal direction perpendicular to its plane. The 
linear response is isotropic in cubic systems but magnetic 
anisotropy develops for non-zero magnetization: the hole 
energy depends on the orientation of M in respect to 
crystal axes and, because of stray field energy Ed, on the 
angle 8 between M and the normal to the film surface. 
A computation of the hole energies E C [M] for relevant 
parameters and [100], [110], and [111] orientations of M 
indicates the that E C [M] can be described by the low- 
est order cubic anisotcopy, so that, taking the stray field 
energy into account E3 

E C (M, 9, if) - E C (M, it/2, 0) = K d (M) cos 2 9 + 
+K cl (M )(sin 4 6 sin 2 tp cos 2 (p + sin 2 9 cos 2 9), (13) 

where K d (M) = 2itM 2 . For 4K d < -K cl the easy axis is 
oriented along [111] or equivalent directions. Otherwise, 
as is usually the case for parameters of (Ga,Mn)As, M 
lies in the (001) plane, and the easy axis is directed along 
[100] for K c i > or along [110] (or equivalent) crystal 
axis in the opposite case. It turns out that the sign of K c \ 
depends on the degree of occupation the hole subbands 
as well as on their mixing by the p-d and kp interactions. 
As a result, the easy axis fluctuates between [100] and 
[110] as a function of p. To quantify the strength of the 
cubic anisotropy in (Ga,Mn)As we computed K c i(M), 
and then the minimum magnitude of an external mag- 
netic field H cu = 2\K cl /M\ (or [i H cu = 2\K cl /M\ in 
the SI units), which aligns the spontaneous magnetiza- 
tion M along the hard direction. Figure ^ shows how 
H cu /M and the direction of the easy axis oscillates as 
a function of p for various Bq = Ap(3M/6giJ,B- As ex- 
pected, H cu tends to zero when Bq decreases. Neverthe- 
less, for \Bq\ — 40 meV (fi M « 85 mT), \i B cu up to 
0.2 T can be expected. Since, however, the orientation 
of the easy axis changes rapidly with p and Bq, intrin- 
sic and extrinsic disorder-which leads to broadening of 
hole subbands-will presumably diminish the actual mag- 
nitude of magnetic anisotropy. 

As shown in Fig. |[ strain has a rather strong influ- 
ence on the valence subbands. It can, therefore, be ex- 
pected that magnetic properties resulting from the hole- 
mediated exchange can be efficiently controlled by strain 
engineering. Indeed, sizable lattice-mismatch driven 
strain is known to exist in semiconductor layers. In 
some cases, particularly if epitaxy occurs at appropri- 
ately low temperatures, such strain can persists even be- 
yond the critical thickness due to relatively high barri- 
ers for the formation of misfit dislocations. We evalu- 
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FIG. 9: Computed minimum magnetic field H cu (divided 
by M) necessary to align magnetization M along the 
hard axis for cubic (unstrained) Gai-^Mn^As film. As a 
function of the hole concentration and the spin-splitting 
parameter Bq, the easy and hard axes fluctuate alterna- 
tively between [110] and [100] (or equivalent) directions 
in the plane of the film. The symbol [110] — > [100] means 
that the easy axis is along [110], so that H cu is applied 
along [100] (Bq — —30 meV corresponds to the satura- 
tion value of M for Gao.95Mno.05As). 



in the region of such low hole concentrations that minor- 
ity spin subbands are depopulated, the easy axis takes 
the [001] direction for the compressive strain, while is in 
the (001) plane for the opposite strain. Such a behav- 
ior of magnetic anisotropy was .recently noted within a 
4x4 model of the valence band.efl However, for the ex- 
perimentally relevant hole concentrations and values of 
Bq (Figs. |l^and 11) the easy axis is oriented along [001] 
direction for the tensile strain, whereas resides in the 
(001) plane for the case of unstrained or compressively 
strained films. This important result, announced already 
in our previous work,tJ is corroborated by the experimen- 
tal study,E3 in which either (Ga,In)As or GaAs substrate 
was employed to impose tensile or compressive strain, re- 
spectively. In particular, for the Gao.o96sMno.o35As film 
on GaAs, for which e xx = —0.2%, the anisotropy field 
HoH un = 0.4 ±0.1 T is observedj£3 in quantitative agree- 
ment with the theoretical results of Fig. O. 

Finally, we mention that strong strain effects may sug- 
gest the importance of magnetostriction in the studied 
compounds. We have not explored this issue yet, and 
note that prior to its examination the question concern- 
ing the collective Jahn- Teller effect in heavily doped p- 
type zinc-blende semiconductors has to be addressed. 



D Optical absorption and magnetic circular 
dichroism 



ate the magnituder-ef resulting effects by using the Bir- 
Pikus hamiltonian,E2l adopted for biaxial strain, as shown 
in Appendix A. Three parameters control strain phe- 
nomena in the valence band: the deformation potential 
b, taken as b =j-^-1.7 eVj£j the ratio of clastic moduli 
C12/C11 = 0. 453^3 and the difference between the lattice 
parameters of the substrate and the layer, Aa. The lat- 
ter is related to relevant components of the strain tensor 
according to, 

£xx = £ yy = Aa/a; (14) 
i-zz = -2e xx ci 2 /cii. (15) 

We have found that biaxial strain has a rather small 
influence on Tq. In the experimentally relevant range 
of hole concentrations 5 x 10 20 > p > 10 20 cm" 3 , both 
tensile and compressive strain diminish Tq- The relative 
effect attains a maximum at p » 2 x 10 20 cm -3 , where 
[T C (e xx) - T C (0)]/T C (P) « -2.4% and -4.9% for e xx = 
1% and —1%, respectively. However, such a strain leads 
to uniaxial anisotropy, whose magnitude can be much 
greater than that resulting from either cubic anisotropy 
or stray fields. The corresponding anisotropy field H un 
assumes the form, 

H un = |2[£ c ([001]) - £ C ([100])]/M + 47rM|[emu];(16) 
H H un = |2[£ c ([001]) - £ C ([100])]/M + Mo M|[SI](17) 

where the last term describes the stray- field effect. Ac- 
cording to results of computations presented in Fig. nG, 



In the case of II- VI diluted magnetic semiconduc- 
tors, detail information on the exchange-induced spin- 
splitting of the bands, and thus on the coupling between 
the effective mass electrons and the localized-spins has 
been obtained from, magnetooptical studies.t3£j Similar 
investigationsOHOS of (Ga,Mn)As led to a number of 
surprises. The most striking was the opposite order of 
the absorption edges corresponding to the two circular 
photon polarizations in (Ga,Mn)As comparing to II- VI 
materials. This behavior of circular magnetic dichroism 
(MCD) suggested the opposite order of the exchange-split 
spin subbands, and thus a different origin of the sp-d in- 
teraction in these two families of DMS. A new light on 
the issue was shed by studies of photoluminescence (PL) 
and its excitatkp..spectra (PLE) in p-type (Cd,Mn)Te 
quantum wells .LtljCj It has been demonstrated that the 
reversal of the order of PLE edges corresponding to the 
two circular polarizations results from the Moss-Burstein 
effect, that is from the shifts of the absorption edges as- 
sociated with the empty portion of the valence subbands 
in the p-type material. This model was subsequently 
applied to interpret qualitatively the magnetoabsorption 
data for metallic (Ga,Mn)AsH3 Surprisingly, however, 
the anomalous sign of the MCD was present also in non- 
metallic (Ga,Mn)As, in athich EPR signal from occupied 
Mn acceptors was seen.E3 Another striking property of 
the MCD is a different temperature dependence of the 
normalized MCD at low-and high photon energies in fer- 
romagnetic (Ga,Mn)AsB3 This observation was taken as 
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FIG. 10: Computed minimum magnetic field H un (di- 
vided by M) necessary to align magnetization M along 
the hard axis for compressive (a) and tensile (b) biax- 
ial strain in Gai—xMn^As film for various values of the 
spin-splitting parameter Bq. The easy axis is along 
[001] direction and in the (001) plane at low and high 
hole concentrations for compressive strain, respectively 

(a) . The opposite behavior is observed for tensile strain 

(b) . The symbol [100] ->■ [001] means that the easy 
axis is along [100], so that H un is applied along [001] 
(Be = —30 mcV corresponds to the saturation value of 
M for Ga .95Mn . 05 As). 
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BIAXIAL STRAIN [%] 

FIG. 11: Computed minimum magnetic field H un neces- 
sary to align the saturation value of magnetization M s 
along the hard axis as a function of biaxial strain com- 
ponent e xx for two values of the hole concentrations in 
Gao.95Mno.05As. The symbol [100] — > [001] means that 
the easy axis is along [100], so that H un is applied along 
[001]. ' 

an evidence for the presence of two. spectrally distinct 
contributions to optical absorptiono 

We begin by noting that according to our two-fluid 
model, the co-existence of strongly and weakly localized 
holes is actually expected on the both sides of the MIT. 
Since the Moss-Burstein effect operates for interband op- 
tical transitions involving weakly localized states, it leads 
to the sign reversal of the MCD, also on the insulating 
side of the MIT. Conversely, the existence of the MCD 
sign reversal can be taken as an experimental evidence 
for the presence of the Fermi liquid-like states on the 
insulating side of the MIT. 

In order to shed some light on those issues we calcu- 
late absorption and MCD spectra in a model that takes 
the complex structure of the valence band into account. 
The band-gap E g is expected to depend on both Mn 
and hole concentration due to the alloy and band nar- 
rowing effects. To take this dependence into account as 
well as to include disorder-induced band tail effectsEj we 
assume, guided by experimental results to be discussed 
below, that E g = 1.2 cV. Our computation of the ab- 
sorption coefficient a(u>) are performed according to a 
scheme outlined in Appendix B, taking the electron ef- 
fective mass m* = 0.07m o , the Kane momentum matrix 
element P = 9.9 x 10~ 8 eVcm, and the refractive in- 
dex n r — 3.5. As shown in Fig. [l^, contributions to a 
originating from particular valence bands are clearly vis- 
ible. Because of the Moss-Burstein shift, the onset and 
the form of a(uj) for particular transitions depend on the 
hole concentration. In particular, a(u>) corresponding to 
the light- hole band exhibits a step-like behavior which, in 
the case of the heavy hole band, is broadened by warp- 
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FIG. 12: Absorption edge in p-GaAs computed for two 
values of the hole concentrations. Inset shows three kinds 
of possible photon-induced transitions corresponding to 
particular valence subbands. Absorption edges associ- 
ated with electron transitions from particular subbands 
are labeled. 

ing. While no quantitative data on a(u>) are available 
for (Ga,Mn)As, the computed magnitude and spectral 
dependence of a(uO reproduce correctly experimental re- 
sults for p-GaAs.O 

The influence of the sp-d band splittings on the ab- 
sorption edge is shown in Fig. |l3|. The computation are 
carried out for the Faraday configuration and with the 
value of the s-d exchange energy aN D = 0.2 eV observed 
in II- VI semiconductors. The theoretical results confirm 
that the Moss-Burstein effect accounts for the sign rever- 
sal of the magnetic circular dichroism (MCD). The en- 
ergy splitting of the absorption edge depends on u but iJts. 
magnitude is similar to that observed experimentally.Ell 
A detail comparison requires from one hand experimental 
information on the absolute values of a(u>) and, on the 
other, more careful consideration of band tailing effects. 
Furthermore, contributions from intra-valence-band and 
from intra-d-shell transitions are expected at low and 
high energy wings of the absorption edge, respectively. 
We predict that not only the former but also the latter 
are substantially enlarged in p-type materials. Indeed, 
the empty valence band states allow for admixtures of 
p-like states to the localized d orbitals. 

The magnetization-induced splitting of the bands is 
seen to lead to a large energy difference between the po- 
sitions of the absorption edges corresponding to the two 
opposite circular polarizations. This may cause an un- 
usual dependence of the low-energy onset of MCD on 
magnetization, and thus on temperature. In particular, 
a standard assumption about the proportionality of MCD 
and magnetization becomes invalid. To find out whether 
the latter is responsible for the anomalous temperature 
dependence of MCD at low photon energies ,E3 we com- 



FIG. 13: Absorption edge for two circular polarizations 
in p-(Ga,Mn)As computed for spin-splitting parameter 
Bq = —10 meV and hole concentration 3.5 x 10 20 cm~ 3 . 
Inset shows how the Fermi see of the holes reverses rela- 
tive positions of the edges corresponding to a + and a~~ 
polarizations. In general, electron transitions from six 
valence subbands contribute to optical absorption. 

pute the differential transmission coefficient that was ex- 
amined experimentally^] 

A = (T+ -T-)/{T + +T~). (18) 

Here T + and T~ are the transmission coefficients for 
right and left circularly polarised light. To take the effect 
of interference into accountj23 these coefficients are cal- 
culated for the actual layout of the samples, which con- 
sisted of a transparent (Ga,Al)As etching stop layer and 
the absorbing (Ga,Mn)As film, each 200 nm thick. The 
same value of the refractive index n r = 3.5 are adopted 
for both compounds. 

Figure |l4| shows the ratio A(w)/A(1.78 eV) computed 
for p — 3.5 x 10 20 cm -3 and various Bq ~ M. In 
the range of high photon energies, u> > 1.6 eV, the re- 
sults collapse into one curve for all values of Bq . This 
means that A(w) = M(T)f(ui) in this range, where 
f(uj) ~ T~ 1 dT/duj. However, in the region of the absorp- 
tion edge, the dependence of A on Bq is by no means 
linear, so that the normalized values A(w)/A(1.78 eV) 
do not follow any single curve. As seen A (a;)/ A( 1.78 
eV) peaks at the greatest value the [Smallest Bq . This is 
the behavior found experimentally.^ We conclude that 
the two observed distinct spectroscopic regionsL3 corre- 
spond to standard band to band transitions, for which 
the proportionality of A to Bq holds, and to the on- 
set of the absorption edge that is shifted and make more 
steep by the Moss-Burstein effect. Actually, the peak val- 
ues of A(w)/A(1.78 eV) determined numerically for the 
low energy region are even greater than that observed in 
(Ga,Mn)As|3 presumably because of scattering broad- 
ening of the absorption edge, neglected in our model. 
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FIG. 14: Spectral dependence of magnetic circular 
dichroism in p-(Ga,Mn)As computed for hole concentra- 
tion 3.5 x 10 20 cm" 3 and various spin-splitting parameters 
Bq. The magnitudes of MCD at given Bq are normal- 
ized by its value at 1.78 eV. 



V. CHEMICAL TRENDS 
A Material parameters 

The ability of the present model to describe success- 
fully various aspects of the ixpjomagnetism in (Ga,Mn)As 
as well as in (Zn,Mn)TeJj'tJE3 has encouraged us to ex- 
tend the theory towards other p-type diluted magnetic 
semiconductors. In this Section, we present material 
parameters that have beep adopted for the computa- 
tions presented elsewhere.B We supplement also the pre- 
vious results^ by the data for (In,Mn)N, (Cd,Mn)S, and 
(Cd,Mn)Se. For concretness, we will assume that 5% of 
the cation sites are occupied by the Mn ions in the 2+ 
charge state, and that the corresponding localized spins 
S = 5/2 are coupled by the indirect exchange interaction 
mediated by 3.5 x 10 20 holes per cm 3 . The enhance- 
ment effect of the exchange interaction among the holes 
in described by the Ferm i liquid parameter Af = 1.2. 
As explained in Sec. II B , no influence of the antiferro- 
magnetic superexchange is taken into account in the case 
of the group III-V and IV semiconductors, in which the 
Mn supplies both localized spins and holes. By contrast, 
it is assumed that in the case of II- VI semiconductors 
x e ff — 0.0297 for x — 0.05 and-JrW? = 1 K, except for 
(Zn,Mn)Te, where Tgp = 2.9 K,Bllf and (Cd,Mn)Te, for 
which T AF = 1.5 KB 

The values of the parameters that are used to deter- 
mined chemical trends are summarized in Tables [| and 
U for cubic and wurzite semiconductors, respectively. If 
no experimentally determined values are available, the 
effective mass parameters are determined by fitting the 
appropriate kp model to results of band structure compu- 
tations. No lattice polaron corrections are taken into con- 



FIG. 15: Energy of p-d exchange interaction for var- 
ious materials containing 5% of Mn as a function of 
lattice parameter. The values shown by solid squares 
were deteraained. Jiom excitonic splittings in the mag- 
netic field ,E3c3'L3'e2I while the empty symbojsjdenote val- 
ues evaluated from photoemission data.olIZl Solid and 
dashed lines represent formulae adopted for the deter- 
mination of the exchange energy for other materials, as 
shown in Tables [I] and 



sideration. Since we are interested in a relatively small 
concentration of magnetic ions, x = 0.05, the effect of 
the Mn incorporation upon the lattice and band struc- 
ture parameters is disregarded. . .. . 

In the case of ZnSe, the recently determinedOlla val- 
ues of the Luttinger parameters ji lead to a negative 
hole mass for the [110] crystallographic direction, an ef- 
fect not supported by the existing theoretical studies of 
the valence band jStructure in this material l 75 li 7s l Accord- 
ingly, an older setE3 of % had been taken for the previous 
calculation.!! The present values of ji are within experi- 
mental uncertainties of the current determinationsO and, 
at the same time, lead to a good description of the-eom- 
puted band structure of the valence band in ZnSe.L3 

In addition to the spin density of states at the Fermi 
level, the Curie temperature is proportional to the square 
of the p-d exchange integral (3. Figure [l5] presents the 
magnitudes of the exchange energy f3N Q as determined 
by photoemission and magnetooptical studies for various 
DMS containing about 5% on Mn. The values of \(3N \ 
are seen to increase-when the lattice parameter decreases. 
This trend stemsLj from the corresponding changes in 
the charge transfer and correlation energies as well as 
from a dependence of the p-d hybridization energy on the 
bond length 6£3 It should, however, be recalled that b, in 
contrast to the average value of the lattice constant, does 
not obey the Vegard low in alloys but rather conserves 
the value corresponding to the end compounds .113 

Thus, in order to obtain the values of j3N a for ma- 
terials, for which no direct determination is available, 
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guided by the results presented in Fig. [n| we assume 
(5N ~ a~ 3 , i.e., (3 — const. More explicitly, for group 
III-V and IV semiconductors we take 

Piii-v = /3(GaMnAs). (19) 

Similarly, for the II- VI materials 

Pu-vi = /3(ZnMnSe), (20) 

where the p-d energy fiNdGaMnAs) — —1.2 eV,0 and 
/37V D (ZnMnSe) = -1.3 eV£3 



B Curie temperatures 

Figure [l6| presents the calculated values of the Curie 
temperature Tc for III-V and II- VI semiconductors con- 
taining 5% of Mn and 3.5 x 10 20 holes per cm 3 . The data 
for Si and Ge are also included. The most remarkable re- 
sult is a strong increase of Tc for materials consisting of 
lighter elements. Actually, Tc exceeding room tempera- 
ture is expected for GaN, InN, and ZnO for the assumed 
values of the Mn and hole concentrations. It has been 
checked for GaN that the value of Tc for the zinc-blend 
modification of this material is by 6% greater than that 
for wurzite case. 

By comparing results of numerical calculations with 
the general Eq. (|7|) for Tc three interrelated mechanisms 
accounting for the chemical trends visible in Fig. |l6| can 
be identified. First, the reduction of the spin density-of- 
states, and thus Tc by the spin-orbit interaction ceases to 
operate in materials with light anions. Second, the effec- 
tive mass, and thus the density-of-states tend to increase 
for materials with stronger bonds. Finally, the smaller 
lattice constant at given x corresponds to the greater 
value of N a x, density of the magnetic ions. It should 
be noted at this point that Tc is proportional to /3 2 , as- 
sumed here to be material independent. This assumption 
corresponds, however, to a strong increase of \8N a \ with 
decreasing lattice constant, as shown in Fig. [l5|. 

It can be expected that the chemical trends established 
here are not altered by the uncertainties in the values of 
the relevant parameters. Our evaluations of the strength 
of the ferromagnetic interactions mediated by the holes is 
qualitatively valid for Mn, as well as for other magnetic 
ions, provided that two conditions are met. First, the 
magnetic electrons stay localized and do not contribute 
directly to the Fermi sphere. Secondly, the holes are de- 
localized and, in particular, do not fopn small magnetic 
polarons, such as Zhang-Rice singletsH 
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FIG. 16: Curie temperatures evaluated for various III- 
V (a) as well group IV and II- VI semiconducting com- 
pounds (b) containing 5% of Mn in 2+ charge state and 
3.5 x 10 20 holes per cm 3 . Material parameters adopted 
for the calculation are displayed in Tables | and [n| 



VI. SUMMARY AND OUTLOOK 



In this paper, theory of ferromagnetism in p-type zinc- 
blende and wurzite semiconductors containing a sizable 
concentration of magnetic ions has been proposed. In has 
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been argued that over the relevant range of the hole den- 
sities the ferromagnetic coupling between the localized 
spins is primarily mediated by delocalized or weakly lo- 
calized holes residing in p-like valence band. Accordingly, 
particular attention has been paid to incorporate into the 
Zener model effects of kp and spin-orbit interactions as 
well as of biaxial strain. It has been demonstrated that 
theory describes qualitatively, and often quantitatively, a 
body of experimental results accumulated over the recent 
years for (Ga,Mn)As. In particular, Curie temperature, 
saturation value of magnetization, hole spin magnetiza- 
tion and polarization, magnetic anisotropies and magne- 
toelastic effects, optical absorption and magnetic circular 
dichroism have been interpreted. 

Giant negative magnetoresistance, sharp field-induced 
insulator-to-metal transition,Eil and a sizable increase pf, 
high frequency conductivity with the magnetic fieldES 
were observed in p-(Hg,Mn)Te. Those findings were at- 
tributed to the growing participation of the light hol es . jn , 
transport when the p-d exchange splitting increases .EiltEl 
This implies a shift of the Drude weight from high to 
low frequencies as a function rfrtte valenG6-.band split- 
ting. Such effects, in both d.c.t-iEa and a.cHo conduc- 
tivity have more recently been detected in (Ga,Mn)As, 
and can qualitatively be interpreted in the same way. 
This provides an additional support for our conclusion 
about the similarity of the mechanisms accounting for 
the hole-mediated exchange interaction in II- VI and III- 
V magnetic semiconductors. However, our work identifies 
also an aspect of ferromagnetism, which points to a dif- 
ference between those two families of magnetic semicon- 
ductors. In the case of II- VI compounds, a short-range 
antifcrromagnetic superexchage lowers the magnitude of 
Tq- This lowering appears to be much less efficient in 
III-V semiconductors, where the Mn ions act as accep- 
tors, compensated partly by donor defects. Thus, the 
localized holes reside preferentially on the Mn pairs, so 
that the hole-mediated ferromagnetic coupling (a vari- 
ant of Zener's double exchange) can overcompensate the 
antifcrromagnetic superexchange. 

The model put forward here has also been used to ex- 
plore the expected chemical trends. It has been found 
that particularly large value of the Curie temperature are 



expected for materials built up from light elements. Im- 
portant issues of solubility limits and self-compensation 
need, however, to be addressed experimentally. In partic- 
ular, the pinning of the Fermi energy by AX- type centers 
or other defects can preclude the increase of the hole con- 
centration in many systems. High pressure research can 
shed some light on this issue. Since, in general, III-V 
compounds can easier be doped by impurities that are 
electrically active, whereas a large quantity of transition 
metals can be incorporated into II- VI materials, a sug- 
gestion has been put forward to,grow magnetic III-V/II- 
VI short period superlattices.E3 Further numerical and 
experimental studies of magnetic semiconductors as well 
as of heavily p-doped non-magnetic systems are expected 
to improve our understanding of the hole-mediated ferro- 
magnetism in zinc-blende and wurzite compounds. This, 
together with exploration of novel quantum structures as 
well as of co-doping and co-alloying, may lead to fabri- 
cation of functional systems. 

On the theoretical side, further work is necessary to 
evaluate quantitative corrections to the mean field theory 
brought about by thermodynamic fluctuations of magne- 
tization. Effects of disorder associated with both random 
distribution of magnetic ions and fluctuations of carrier 
density near the metal-insulator transition are other open 
issues. In particular, unknown is the nature of evolu- 
tion of static and dynamic magnetic phenomena on ap- 
proaching the strongly localized regime. The above list 
of interesting problems is by no means exhausting. With 
no doubt we will soon witness many unforeseen develop- 
ments in the field of carrier-mediated ferromagnetism in 
semiconductors . 
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APPENDIX A: EFFECTIVE-MASS HAMILTONIANS 

The purpose of this Appendix is to, provide the explicit form of the effective mass hamiltonian that, in addition 
to the standard kp and strain terms,E3 contains a contribution of the p-d exchange interaction in the molecular field 
approximation. The latter constitutes a generalization of the previous approacheafjCj by allowing for the arbitrary 
orientation of the magnetization M with respect to the crystal axes. 

Zinc-blende semiconductors are considered first. We take explicitly into account four Tg and two T7 valence 
subbands, for which we chose the basis functions in the form: 

ui = —7=(X + iY) I, (Al) 
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U2 = i-L[(X + iF)|-2ZT], (A2) 

u 3 = -j=[(X-iY)1+2Zl], (A3) 

U4 = i _L(X-iF)|, (A4) 

U5 = -L[(X + iF)|+ZT], (A5) 

U6 = i4=[-(X-iy)T+^|], (A6) 



where X, Y, and Z denote Kohn-Luttinger amplitudes, which for the symmetry operations of the crystal point group 
transform as p x ,p y , and p z wave functions of the hydrogen atom. 

In the above basis the corresponding Luttinger-Kohn matrices assume the form: 



1 kp matrix 



Hkp — —- 





' P + Q 


L 


M 





iL/y/2 


-\\f2M 




L* 


P-Q 





M 


-iy/2Q 


iv4V2L 




M* 





P-Q 


-L 


-\y/ij2L* 


-iV2Q 


2m 





M* 


-L* 


P + Q 


-iM* 


-iL*/V2 




-LL*/\/2 


iV2Q 


iy/3/2L 


iv^M 


P + A 







iV2M* 




iV2Q 


iL/y/2 





P + A 



Here, 



P = ^k 2 , 

Q = 72(k 2 x + k 2 y -2k% 
L = -i2V3j3(k x - ik y )k z , 
M = V3[j2(k 2 x -k 2 y )-i2 l3 k x k y ], 
A = 2m A /h 2 . 



(A7) 



(A8) 
(A9) 

(A10) 
(AH) 
(A12) 



2 p-d exchange matrix 



H P d = B G 



Here, 



3b x w z i^3b x w- 
-i\/3b x w + (2b z - b x )w z 

—2\b z w + 


\/6b x w + -iV2(b x + b z )w z 





2ib z w- 

-(2b z -b x )w z \\fZb x w- 

—i\/3b x w+ —3b x w z 

V2b z w- 







V2b z w + iV2(b x + b z )w z -^b x w. 



B G =l3M/6gp B , 

w z = M z /M, 

w± = (M x ± iM y )/M, 

b z = x /0, 

b x = x /0, 



V6b x W- 
\\f2(b x + b z )w z 
V2b z w + 


-(26a; - b z )w z 
-\b z w + 





— \[2b z W- 
-\\f2{b x + b z )w z 
-\/&b x w + 
ib z w- 
(2b x - b z )w z 



(A13) 



(A14) 
(A15) 
(A16) 
(A17) 
(A18) 



where in cubic materials b z = b x = 1. 
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3 Biaxial strain matrix 



H bs = b 



Here, b is the deformation potential and 



" -Qe 





R ( 








-iv^i?. 





Q, 







iV2Q e 





Re 





Q, 













Re 





-Qe 


-iV2R e 








-iV2Qe 





lVZRe 













-iV2Qe 












Re 



V3(e xx - €yy)/2. 



(A19) 



(A20) 
(A21) 



Since we are interested in the effect of the biaxial strain in the (001) plane, only the terms involving the diagonal 
components en of the deformation tensor are included. For the same reason, we allow for the corresponding anisotropy 
of the exchange integrals (3 X = (3 y ^ f3 z , though we expect that to a good accuracy f3 x — [3 Z in real systems. The 
latter is assumed in the main body of the paper. 

In the presence of the magnetic field B the Luttinger-Kohn kp matrix is a sum of the Zeeman and Landau-parts, 
Hk P = Hz + Hl, where Hl is written below for B || [001] and neglecting some terms proportional to 72 — 73O 



H, 



3kw z /2 
— iv3fttf+ /2 





i^/Snw- /2 

KW z /2 







lKW- 
— KW z /2 

— i^/3kw+ /2 



-VQ(k + 1)w+/4 iV2(k+1)w z /2 -V2(k + l)io_/4 
V2(k + 1)w + /4 -iV2(k + 1)w z /2 





_ 

iy3nW-/2 
—3kw z /2 


V6{k + l)w-/4 



-\/6(« + l)tu-/4 
-iy/2(K+l)w z /2 V2(k + 1)w_/4 
-V2(k + 1)w+/4 iV2(n+l)w z /2 





(k+1/2)w z 
i(n + l/2)w+ 



-\(k + 1/2)w- 
-(k+1/2)w z 



where g Q is the Lande factor of the free electron, and 



w z 
w± 



B z /B 
B T ± iB n , 



For the Landau quantum number n > 



Hr 



2m 



hh + s(7i +72)(2n- l) b^/n 
-by/nh + s(ti - 72) (2n +1) 
Cv/n(n + 1) 

cV(" + l)(n + 2) 

ib^/n/2 iq + iy/2j 2 s{2n + 1) 

ife v /3(n+l)/2 



cy n(n + 1) 
cVTn + l)(n + 2) 
+ s(7i -7 2 )( 2n + 3) 



ic v / 2(n+T)n iOV3(,n + lj/2 ig • 


-ig-iV2s7 2 (2n+l) 

-by/~(n + 2) 

h h + (71 + 7a)g(2n + 5 ) -ic^2(n + l)(ra + 2) 

ic- v /2( n + 2)(n+ 1)) /i s + 7 X s(2n + 1) 

i6V(" + 2)/2 



i&v/3(n+ l)/2 
■iy27 2 s(2n + 3) 



i& A /n/2 
i6y/ 3(n+ l)/2 
iby/3(n+ l)/2 



-ic-y/2n(n + 1) 

-ig - i V27 2 s(2n + 3) 
iV(^ + 2)/2 


/i s +7is(2n + 3) 



(A22) 



(A23) 
(A24) 



(A25) 
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Here. 



s = 



eB/h, 
h h = ( 7l -2 l2 )k 2 z , 
hi = ( 7l +2 l2 )k 2 z , 

-2V2 72 fe z 2 , 
-i2\/6s 7 3fc z ; 
V / 3s( 72 + 73 )- 



/la = 

q : 

b 

c = 



If n = 0, the wave function basis does not contain the Ui term, so that 



2m 



v /(n + 2)(n+l) 



/i; + s( 7 i - 72 )(2n + l) 

/i; +s( 7 i - 72 )(2n + 3) 

h h + ( 7 i + 72 )-s(2n + 5) 
i Cy /2( n + 2)(n+ 1) 
iV(n + 2)/2 



cV(^ + 2)(n+ 1) 
ig + i V2 72 s(2n + 1) 
i6^3(n+l)/2 



n - 72 j 

by/n + 2 
i&V3(n + l)/2 



ig + iV2 72 s(2n + 3) 

-ig - i V2g 72 (2n + 1) 16^/3(71 + l)/2 



ib ^/3(n+ l)/2 -ig - i V2 72 s(2n + 3) 
-ic v /2(n + 2)(ra+ 1) i6^/(n + 2)/2 
/i s + 7i s(2n + 1) 

/i s + 7 is(2n + 3) 

If n = — 1, the wave function basis does not contain the u\, u 2 , and u 5 terms, so 

h 2 \ hi + s( 7 i - 72 )( 2n + 3) -b^Jn + 2 -ig - K/2 72 s(2n + 3) 



' that 



H L = - 



2m 



hi + s( 7 i - 72 )( 2n + 3) -b^n + 2 

bVnT2 h h + ( 7 i + 72 )g(2n + 5) \byj(n + 2)/2 

+ iV2 72 s(2n + 3) i6^(n + 2)/2 /i s + 7 is(2n + 3) 



lg- 



If n=-2, the wave function basis contains only 114 term, so that 



H L = 



2m 



-[hh + (li + 72 )s]. 



(A26) 
(A27) 
(A28) 
(A29) 
(A30) 
(A31) 
(A32) 



(A33) 



(A34) 



(A35) 



the 



The solution of the eigenvalue and eigenfunction problem for the hamiltonian H = H^ p + H, p( i + Hi, s or, in 

magnetic field, H = Hl + Hz + H p d + Hi, s gives the hole energies £j and eigenvectors for the six relevant hole 
subbands. 



4 Wurzite compounds 

For wurzite compounds, we chose the basis in the form, 





(A36) 


Ui = ^=(X-iY)U 


(A37) 
(A38) 
(A39) 


u 5 = -^=(X + iY)l 

"6 = Z i, 


(A40) 
(A41) 
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In the above basis, the kp matrix reads 



2m 



Here. 



F 


-K* 




-H* 








-K 


F-2A 2 




H 





V2A3 


-H 


H* 


L 


- Ai - A 2 


V^Ag 













F 


-A" 


H 


n 
U 


n 
u 




V <ZL±3 "A 


_r — ZZA2 


rr* 
— / / 


o 


./OA, 




77"* 

w 11 


—H L 


- Ai - 






L 


= A 1 k 2 z +A 2 {kl -l 










T 


= A 3 k 2 z +Ai(kl -I 










F 


= L + T, 










k + 


A^a; + iky, 










K 


= A 5 k 2 + , 










H 


= A 6 k + k z , 







where Aj and -Aj are the valence band splittings and kp parameters, respectively, 
p-d exchange matrix is given by 



H P d = 3Bq 



b x w z 











Mm 


" 





b x w z 





b x w m 














b z w z 








b z w m 





b x w p 





~b x w z 








b x w p 











-b x w z 











b z Wp 








-b z w z _ 



(A42) 



(A43) 
(A44) 
(A45) 
(A46) 
(A47) 
(A48) 



(A49) 



APPENDIX B: NUMERICAL PROCEDURE 



We aim at evaluating the Helmholtz free energy F c at given hole concentration p as a function of the Mn magneti- 
zation M, 



F c {p, M) = -k B T J &eN{e) In (1 + cxp[-(e(M) - e F (p, M))/k B T]) + pe F (p, M), 



(Bl) 



where the thermodynamic density-of-states N(e) — dp/de. If T < 200 K then in the studied cases the hole liquid is 
degenerate, so that F c assumes a simple form, 

rp 

(B2) 



F c (p,M)= / d P /e(M,p/). 
Jo 

Thus, in order to obtain F c (p,M) we have to integrate the dependence of the partition function or of the Fermi 
energy on the hole concentration for various M. The actual calculation proceeds in a standard way. First, the length 
of the wave vectors ki{6, if) for each of six hole subbands is determined by solving the inverse eigenvalue problem for 
given values of e, polar angle 9 and azimuth angle ip in the fc- vector space. Then, the hole concentration is obtained 
by an integration of J^i kf(6, (p)/24ir 3 over cos (9 and ip. 

In order to calculate the contribution M c of the hole magnetic moments to the total magnetization, the eigenvalue 
problem is solved directly making it possible to determine the Gibbs thermodynamic potential G c , 



G C = ~Y^ k B T(eB/2nh)\n{l + exp[-{e l {n,k z ,M) - e F )/k B T}) . 



(B3) 



We calculate the absorption coefficient a for the two circular polarizations cr^ , taking into account k conserving 
electron transitions from the six valence band subbands, index i, to the two spin branches of the conduction band, 
index j, assuming that the hole liquid is strongly degenerate. In such a model J^Do 



, 47r 2 e 2 P 2 v-^ 
Q =tz 2^ 



r, 2 



cn r u> 



1,3 



dcos6> 
1 Jo 



2 " k %3 {uj)m { r 3) {uj) , 2 
dip 8^ |M ^' ■ 



(B4) 
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The wave vectors kij(u),9,ip) corresponding to twelve possible transitions are determined by the energy conservation 
and the position of the Fermi level. The joint density of states effective mass corresponding to these fcjj is given by 



TO fe) = (l/m e - 



1 d£j 



(B5) 



The matrix elements M^(6, ip) for the two light polarizations and involving electron transitions to the spin down and 
spin up conduction subbands are given by 

(B6) 
(B7) 
(B8) 
(B9) 





(4) 
= a i > 






= afVv / 3 + iaf\ 




Ma 


= a^/V3-ia!f\ 




M~ 2 







where a" is the n-th component of the eigenvector corresponding i-th valence subband at ip). 
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APPENDIX C: MATERIAL PARAMETERS 



In Tables I and II we summarize band structure pa- 
rameters of parent compounds employed for the evalu- 
ation of chemical trends in cubic and wurzite magnetic 
semiconductors . 



a the value determined from Eq. ( [Til ) 

h the value determined from the data for the wurzite structure 
(Table 0) 

c K. Kim, W.L.R. Lambrecht, B. Segall, and M. van 
Schilfgaarde, Phys. Rev. B 56, 7363 (1997); in this paper 
the Dresselhaus parameters L,M,N are denoted as A,B,C, 
respectively; see also, I. Stolpe, N. Puhlmann, H.-U. Miiller, 
O. Portugall, M. von Ortenberg, D. Schikora, D.J. As, B. 
Schottker, and R. Lischka, Physiea B 256-258, 659 (1998) 

the value from photoemissionEfl 
e the value determined from Eq. ( po| ) 
^ see the main body of the text 

9 H.W. Holscher, A. Nothe, and Ch. Uihlein, Phys. Rev. B 
31, 2379 (1985) 
h magnetoreflectivit\ 
% magnetoreflectivity 
3 magnetoreflectivity 



a Calculated from Eq. (J19J) 
b Calculated from Eq. ( pCp ■— ■ 
c interpretation of rrmgnetooptical dataE3 
d magnetoabsorptionu 



TABLE I: Material parameters of selected cubic semi- 
conductors and the values of the p-d exchange energy 
j3N . Except for the parameters for which references are 
provided, the values of the lattice constant a a , spin-orbit 



splitting A , and . 
Landolt-BdrsteinX 



parameters 7$ are taken from 





o D (A) A (eV) 


7i 


72 


73 


PN (eV) 


Si 


5.43 


0.044 


4.285 


0.339 


1.446 


-1.35 a 


Ge 


5.66 


0.29 


13.38 


4.24 


5.69 


-1.20 a 


A1P 


5.47 


0.1" 


3.47 


0.06 


1.15 


-1.33 a 


AlAs 


5.66 


0.275 


3.25 


0.64 


1.21 


-1.19 a 


GaN 


4.50 b 


0.018 6 


2.463 c 


0.647 c 


0.975 c 


-2.37 a 


GaP 


5.45 


0.080 


4.05 


0.49 


1.25 


-1.34 a 


GaAs 


5.65 


0.34 


6.85 


2.1 


2.9 


-1.2 d 


GaSb 


6.09 


0.76 


13.3 


4.4 


5.7 


-0.96° 


InP 


5.87 


0.108 


5.15 


0.94 


1.62 


-1.07 a 


InAs 


6.06 


0.38 


20.4 


8.3 


9.1 


-0.98 a 


ZnS 


5.401 


0.070 


1.77 a 


0.30 a 


0.62 a 


-1.5 e 


ZnSe 


5.67 


0.43 


2.95 ; 


0.6 / 


1.11 s 


-1.3 h 


ZnTe 


6.10 


0.91 


3.8 


0.72 


1.3 


-i.r 


CdTe 


6.48 


0.95 


4.14 


1.09 


1.62 J 


-0.88 J 



TABLE II: Material parameters of selected wurzite semi- 
conductors. The values of the lattice parameter a a = 
(v / 3a 2 c) 1 / 3 and the splittings Aj of the valence band at 
the r poinlJii II- VI compounds are taken from Landolt- 
Borsteini^sEB Theoretical studies provide the parameters 
Ai &iuL.A l for GaNjEil InNjE3 and A l for ZnO, CdS, and 
CdSeO 

Material GaN InN ZnO CdS CdSe 



a (A) 
Ai (eV) 
(eV) 
(eV) 



At 
M 
A 3 
Ai 
A 5 
As 

(3N (eV) 



4.503 4.974 
0.0036 0.017 
0.005 
0.0059 
-6.4 



-0.5 

5.9 
-2.55 
-2.56 
-3.06 
-2.37 a 



0.001 
0.001 
-9.28 
-0.6 
8.68 
-4.34 
-4.32 
-6.08 
-1.76 a 



4.567 
0.04 
0.0 
0.0 
-2.41 
-.044 
2.11 
-1.06 
-1.06 
-1.51 
-2.48 6 



5.845 
0.030 
0.022 
0.022 
-5.92 
-0.70 
5.37 
-1.82 
-1.82 
-1.36 
-1.55 c 



6.078 

0.039 

0.139 

0.139 

-10.2 

-0.76 

9.53 

-3.2 

-3.2 

-2.31 

-1.27 d 
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